Constraints on effective constitutive parameters of certain 
bianisotropic laminated composite materials 



AKHLESH LAKHTAKIA 

Nanoengineered Metamaterials Group 
Department of Engineering Science & Mechanics 
212 Earth & Engineering Sciences Building 
Pennsylvania State University 
University Park, PA 16802-6812 
USA 

Abstract 

When the electrically thin unit cell of a laminated composite material is made of two 
bianisotropic sheets whose constitutive properties in the thickness direction are decoupled 
from the constitutive properties in the interfacial planes, the laminated composite material 
can be homogenized into a material not all of whose constitutive parameters are independent 
of each other. This non-independence of the constitutive dyadics of the constituent materials 
and the homogenized composite material is captured by two simple constraints, which may not 
hold if even one of the two constituent materials has more complicated constitutive properties 
than stated above. 
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1 Introduction 

The electromagnetic response properties of laminated composite materials have been and continue 
to be of technoscientific importance (Herpin, 1947, Abeles, 1950; Lafait et at, 1990; Neelakanta, 
1995; Lakhtakia, 1996). A laminated composite material is made by stacking together sheets of 
different materials one on top of another, in order to form a unidirectionally nonhomogeneous 
material with piecewise-uniform constitutive properties. When a monochromatic plane wave is 
incident on this stratified composite material, standing waves are set up in each sheet, and elec- 
tromagnetic energy may also be absorbed therein. In addition, the phenomenons of reflection and 
transmission occur. Across any interface, the tangential components of the electric and the mag- 
netic field phasors are continuous. These tangential components are the essential ingredients of a 
two-point boundary value problem that can be adequately stated and solved using matrix calculus. 
Whereas 2x2 matrixes suffice for isotropic dielectric-magnetic sheets, 4x4 matrixes are needed for 
sheets with more complicated electromagnetic properties (Teitler & Henvis, 1970; Mrozowski, 1986; 
Lakhtakia, 1987). 

A homogenization problem emerges when electrically thin sheets are stacked periodically and 
the unit cell is also electrically thin (Wiener, 1912; Rytov, 1956; Rumsey, 1964). Such a period- 
ically nonhomogeneous composite material whose unit cell comprises two or more sheets may be 
considered equivalent to a homogeneous material that necessarily has direction-dependent consti- 
tutive properties. Analytically straightforward techniques can be employed to predict the effective 
(i.e., post-homogenization) constitutive parameters of a laminated composite material (Reese & 
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Lakhtakia, 1991; Ramakrishna & Lakhtakia, 2009). The effective constitutive parameters depend 
on the constitutive parameters of the constituent materials as well as on their volume fractions. 

Are all effective constitutive parameters necessarily independent of each other? Or, knowledge 
of some effective constitutive parameters is sufficient to determine the others without recourse to 
the volume fractions of the constituent materials. This brief communication addresses this issue. 

An cxp(— iut) time-dependence is implicit, with uo denoting the angular frequency. The perme- 
ability and permittivity of free space (i.e., vacuum) are denoted by /i and s a , respectively. Vectors 
are in boldface, dyadics are underlined twice, column vectors are in boldface and enclosed within 
square brackets, and matrixes are underlined twice and similarly bracketed. Cartesian unit vectors 
are identified as u x , u y and u z , with the z axis oriented normal to the sheets. 



2 Analysis 



Let us consider a periodic laminated composite material whose unit cell is made of two sheets 
of dissimilar, homogeneous, bianisotropic materials labeled 1 and 2. The Tellcgcn constitutive 
relations of the two materials are written in the frequency domain as (Lakhtakia, 1987) 



D = e a (|(") • E + • H) 
B = n (>>) -H + £ (n) -e) 



ne{l,2} . 



(1) 



The relative permittivity dyadic and the relative permeability dyadic of the n-th material are de- 
noted by e(") and fjM 1 ^ , respectively, whereas the dyadics and (3^ denote the magnetoelectric 
properties. The thickness of the sheets made of the n-th constituent material is denoted by d n , so 
that 

fn = d n /(d 1 +d 2 ), ne{l,2}, (2) 

is the volume fraction of the n-th material, with fi = 1 — f\ . 

Wave propagation in the laminated composite material can be handled by using the spatial 
Fourier transform as (Krowne, 1984) 



E(x,y,z) — J„ 27r e(z, k, ip) exp [iK,(xcosip + y simp)] dipdn 



H(x,y,z) = h(z, n,ip) exp [in(x cos ip + ysinip)] dip dn 



(3) 



where n is the spatial frequency and ip is an angle. Substitution of this representation along with 
the constitutive relations (1) in the two Maxwell curl equations leads to the 4x4-matrix ordinary 
differential equation (Lakhtakia, 1987) 



^[F(z, K ,iP)}=i[E {n) ( K ^) 



\F(z,k,i/>)] , 

in any sheet made of the n-th material. The column vector 



e{i,2} 



(4) 



[F(z,k,iP)] 



e x (z,n,ip) 
e y (z,K,ip) 
h x (z,K,ip) 

hy(Z 7 K, Ip) 



(5) 



represents components of the electric and magnetic field phasors that are tangential to the bi- 



material interfaces in the laminated composite material. The 4x4 matrix 



£ (n W) 



is too 
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cumbersome to reproduce here, but it can be put in the following form (Lakhtakia & Weiglhofer, 
1997): 



P [n) M) 



(n) 



+ 



£ 3 



(n) 



;^}+« 2 Ig^W] , nG{l,2}. (6) 



Provided the sheets are electrically thin (Lakhtakia & Krowne, 2003; Mackay, 2008), the lami- 
nated composite material is equivalent to a homogeneous material with the following constitutive 
relations: 

D = e a (g e ff ■ E + g?ff ■ H) 

B //..(//'-' Ml • i -' -K) ' ' 



Equations (3) and (4) govern wave propagation in the homogenized composite material, except 



that the matrix 



P {n) (K, V) 



must be replaced by 



eff 



+ 



eff 



eff 



=4 



(8) 



As the sheets and the unit cell are electrically thin, invocation of the long-wavelength approxi- 
mation (Reese & Lakhtakia, 1991; Lakhtakia & Krowne, 2003) yields the relation 



£ e// ( K ,V)] = h [z (1) (^)] +/ 2 [z^^vo] , V{ K , V} 



(9) 



which can be used to determine |e e -^, [f^, a e ^ ', /3 e ^ j in terms of |e^\ fj,^\ j, 
{| (2) , M (2) , ^ (2) , P (2) }, fu and f 2 . The equation 



is satisfied \A/> £ [0, 2tt) by 



Q e / f W] =/i fe^Wl +/ 2 fe 2) (V>) 



e ^// =7 (/ 1 ei 1 J<I>(^ 2 )+/ 2e i 2) <i>( 1 - 1 )) 
^ / =7(/i^ ) <f (2 ' 2) +.f 2 4 2 , ) * (M) ' 

^ / =7(/iM 1 ^ (2 ' 2) +/2 4 2 ^ (1 ' 1) 



(10) 



(11) 



where = u z ' 1 u z, etc., and 



7 -l = f 2 $(2,2) + /2 $(1,1) + /i/2 ( $( 1,2) + $(2,1)) 

Solution of the equations 



Q 



eff 



h Q {1) +h Q 



)(2) 



G {1,2,3}, 



(12) 



(13) 



yields analytical expressions for the remaining components of |e e ^ i , fj, e f f , a e ^ f , (3 e ^ j , but those 
expressions are far too unwieldy, in general, for reproduction here. 
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2.1 A special case 



A special case emerges when all four constitutive dyadics of both constituent materials are of the 
form 

i (n) 

u z • A (n) = u, = 



so that a Z x = 
simplifications 



An) 



(n) 



,(«) _ 



'■yz 



= 0. Then, 



Q 



(n) 



(n) 



= a % f = <fJJ = <JJ = i a e {e, /z, a, /?} 
follow from eqs. (13). Furthermore, those equations yield 



ne{l,2}, (14) 

: [£],("€ {1,2}), and the 
(15) 



,e// 



ha$+ha$, a e {e,/j,,a,(3} , pe{x,y}, q e {x,y} . 



,(2) 



(16) 



These expressions hold in addition to eqs. (11) for a e JJ . 

Equations (11) and (16) indicate that — once e e J q ^ and [fj q ' f , (p G {x,y, z} , q E {x,y, z}), have 
been calculated by using the constitutive parameters and the volume fractions of both constituent 
materials — a e J^ and ffi^f can be obtained without using the volume fractions at all. Thus, 10 of 
the 20 effective constitutive parameters in the special case under consideration require knowledge 
of the volume fractions, but thereafter the remaining 10 effective constitutive parameters do not. 
Accordingly, all effective constitutive parameters are not independent of each other. 



2.2 General consideration 

Equations (11) and (16) suggest the formulation of the constraints 



Dct 



r(l) g(2) 

(1) „.(2) 



a 



eff 



= Det 



r(l) g(2) 

,(1) „.(2) 



re// 
eff 



lF> ~P {2) li eff 



= 0. 



(17) 



Both certainly hold true for the special case treated in Sec. 2.1. Could these constraints hold in a 
more general sense? 

When even one of the two constituent materials is more general than in Sec. 2.1, analytical 

expressions for je 6 -^, , a e ^ , /3 e ^ j turn out to be so huge that analytical expressions of the 

determinants in the two constraints (17) could not be manipulated to ascertain if both constraints 
hold true in general. 

Several numerical experiments were conducted, wherein all 36 effective constitutive parameters 
were computed with fi G (0, 1) and with none of the 72 constitutive parameters of the two con- 
stituent materials of null value. The determinants in eqs. (17) turned out be significantly different 
from zero, indicating thereby that the constraints (17) do not hold in general. 



2.3 Unit cell with 3 or more sheets 

Suppose the unit cell of a certain laminated composite material comprises three different sheets. 
Two adjacent sheets can be homogenized into one thicker sheet, and the conclusions in Sec. 2.1 and 
2.2 apply to that homogenization. This thicker sheet and the third sheet can also be homogenized 
into a single sheet of the thickness of the unit cell, and the conclusions in Sec. 2.1 and 2.2 also 
apply to that homogenization. This procedure can be adopted for unit cells with 4 or more sheets, 
so long as the unit cell, after all homogenization steps, remains electrically thin. 
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3 Concluding Remarks 



Suppose that the unit cell of a laminated composite material is made of two bianisotropic sheets 
whose constitutive properties in the thickness direction arc decoupled from the constitutive prop- 
erties in the intcrfacial planes. Provided that the long-wavelength approximation is applicable, 
the laminated composite material can be homogenized into a material all of whose constitutive 
parameters arc not independent of each other. This non-independence of the constitutive dyadics 
of the constituent materials and the homogenized composite material is captured by two simple 
constraints. When even one of the two constituent materials has more complicated constitutive 
properties, the constraints are not expected to hold. 

Given that magnetoelectric properties are considerably rarer and usually weaker than dielectric- 
magnetic properties, the two constraints, when valid, should be considered as applicable on the 
magnetoelectric dyadics. Isotropic chiral and biisotropic materials, general uniaxial bianisotropic 
materials, gyrotropic materials such as ferrites and Faraday chiral materials, etc., (Mackay & 
Lakhtakia, 2008) all lie within the scope of the two constraints. In the realm of classical electro- 
magnetics (i.e., at frequencies not exceeding about 750 THz), we can expect the homogenization 
procedure to hold for sheets as thin as 10 molecular diameters (Kim et at, 2005). Thus, the con- 
clusions drawn in this communication can also be expected to hold at optical and lower frequencies 
for many types of nanotextured thin films. 

Acknowledgments. Thanks are due to S. Anantha Ramakrishna (USER, Mohali, India) for 
discussions. The Charles Godfrey Binder Endowment at Penn State is gratefully acknowledged for 
partial financial support. 
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